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The basic notion of percolation in physics assumes the emergence of a giant 
connected (percolation) cluster in a large disordered system when the den- 
sity of connections exceeds some critical value. Until recently, the percolation 
phase transitions were believed to be continuous, however, in 2009, a remark- 
ably different, discontinuous phase transition was reported in a new so-called 
"explosive percolation" problem. Each link in this problem is established by 
a specific optimization process. Here, employing strict analytical arguments 
and numerical calculations, we find that in fact the "explosive percolation" 
transition is continuous though with an uniquely small critical exponent of the 
percolation cluster size. These transitions provide a new class of critical phe- 
nomena in irreversible systems and processes. 
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The modern understanding of disordered systems in statistical and condensed matter physics 
is essentially based on the notion of percolation (1). When one increases progressively the 
number of connections between nodes in a network, above some critical number (percolation 
threshold) a giant connected (percolation) cluster emerges in addition to finite clusters. This 
percolation cluster contains a finite fraction of nodes and links in a network. The percolation 
transition was widely believed to be a typical continuous phase transition for various networks 
architectures and space dimensionalities (2), so it shows standard scaling features, including a 
power-law size distribution of finite cluster sizes at the percolation threshold. Recently, how- 
ever, it was reported that a remarkable percolation problem exists in which the percolation 
cluster emerges discontinuously and already contains a finite fraction of nodes at the percola- 
tion threshold (3). This conclusion was based on numerical simulations of a model in which 
each new connection is made in the following way: choose at random two links that could be 
added to the network, but add only one of them, namely the link connecting two clusters with 
the smallest product of their sizes. To emphasize this surprising discontinuity, this kind of per- 
colation was named "explosive" (3). Further investigations of "explosive percolation" in this 
and similar systems, also mainly based on numerical simulations, supported this strong result 
but, in addition, surprisingly for abrupt transitions, revealed power-law critical distributions of 
cluster sizes (4-10) resembling those found in continuous percolation transitions. This self- 
contradicting combination of discontinuity and scaling have made explosive percolation one of 
the challenging and urgent issues in the physics of disordered systems. 

Here we resolve this confusion. We show that there is not actually any discontinuity at the 
"explosive percolation" threshold, contrary to the conclusions of the previous investigators. We 
consider a simple representative model demonstrating this new kind of percolation and show 
that the "explosive percolation" transition is a continuous, second-order phase transition but, 
importantly, with an uniquely small critical exponent j3 0.0555 of the percolation cluster 
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size. 

One of the simplest systems in which classical percolation takes place is as follows. Start 
with iV unconnected nodes, where N is large, and at each step add connection between two 
uniformly randomly chosen nodes. In essence, this is a simple aggregation process (11), in 
which at each step, a pair of clusters, to which these nodes belong, merge together (Fig. 1A). If 
we introduce "time" t as the ratio of the total number of added links in this system, L, and its 
size N, i.e., t = L/N, then the percolation cluster of relative size S emerges at the percolation 
threshold t c — 1/2 and grows with t in the following way: S ~ 5 13 , where 5 = \t—t c \ and (3 = 1. 
At the critical point, the cluster size distribution (fraction of finite connected components of s 
nodes), n(s), in this classical problem is power-law: n(s) ~ s~ T with exponent r = 5/2, see 
Ref. (1). 

In this report, we consider a direct generalization of this process. Namely, at each step, we 
sample twice: 

(i) choose two nodes uniformly at random and compare the clusters to which these nodes 
belong; select that node of the two ones, which belongs to the smallest cluster; 

(ii) choose a second pair of nodes and, again, as in (i), select the node belonging to the 
smallest of the two clusters; 

(iii) add a link between the two selected nodes thus merging the two smallest clusters (Fig. IB). 

Repeat this procedure again and again. Note that in (i) and (ii) a cluster can be selected several 
times. This is the case for the percolation cluster. These rules contain the key element of other 
explosive percolation models, e.g., model (3), namely, for merging, select the minimal clusters 
from a few possibilities. Importantly, our procedure provides even more stringent selection of 
small components for merging than model (3) since guarantees that the product of the sizes of 
two merging clusters is the smallest of the four possibilities (each of the first pair of chosen 
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nodes (i) may connect with any node of the second pair (ii)) in contrast to selection from only 
two possibilities in model (3). Consequently, if we show that the transition in our model is 
continuous, than model (3) also must have a continuous transition. More generally, in each 
sampling, one can select at random m nodes thus choosing the minimal cluster from m clusters. 
Classical percolation corresponds torn = 1 (Fig. 1A). In this report we mostly focus on m = 2 
(Fig. IB). One should stress that the explosive percolation processes are irreversible in stark 
contrast to ordinary percolation. In the latter, one can reach any state either adding or removing 
connections. For explosive percolation, only adding links makes sense, and an inverse process 
is impossible. 

We simulated this irreversible aggregation process for a large system of 2 x 10 9 nodes. When 
plotted over the full time range, the obtained dependence S(t) shows what seems to be disconti- 
nuity at the critical point t c (Fig. 1C) similar to previous results, but a more thorough inspection 
of the critical region (log-log plot in Fig. ID) reveals that the obtained data is definitely better 
fitted by the law aS 13 , which indicates a continuous transition, than, say, by 0.3 + b5@. Fitting 
the data of our simulation by the law So + b5@, we find that S is at least smaller than 0.05. This 
shows that for a definite conclusion, even so large a system turns out to be not sufficient, and 
a discontinuity can be ruled out or validated only by analytical arguments for the infinite size 
limit. 

We address this problem analytically and numerically by considering the evolution of the 
size distribution P(s) for a finite cluster of s nodes to which a randomly chosen node belongs: 
P(s) = sn(s)/{s), where (s) is the average cluster size (the ratio of the number of nodes and 
the total number of clusters). This distribution satisfies the sum rule J2 S P( s ) = 1 — 5'. Another 
important characteristic in this model is the probability Q(s) that if we choose at random two 
nodes than the smallest of the two clusters to which these nodes belong is of size s. Q(s) 
provides us with the size distribution of merging clusters. Here J2 S Q( S ) — 1 — S 2 . If we 
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introduce the cumulative distributions P cum (s) = J2^= S P( U ) an d Qcum = J2T=sQ( u )> men 
probability theory gives Q cum (s) + S 2 = [P cum (s) + S] 2 . Consequently 

Q(S) = [Pcum(s) + Pcum(s + 1) + 2S]P( S ) = [2 - 2P(l) - . . . - 2P( S - 1) - P(s)]P( S ), (1) 

that is Q(s) is determined by P(s') with s' < s. The evolution of these distributions in the 
infinite system is exactly described by the master equation: 

^M = S £ Q(«,f)QM)-2sQ(M), (2) 

which generalizes the standard Smoluchowski equation. We derived Eq. (2) only assuming 
the infinite system size. The only difference from the classical percolation problem (11) is the 
presence of the distribution Q(s,t) instead of P(s,t) on the right-hand side of this equation. 
Thus we have a chain of coupled equations, which should be solved analytically or numerically. 
To find a numerical solution, first solve the first equation of the chain, which gives P(l,t). 
Substitute this result into the second equation and solve it, which gives P(2, t), and so on. In this 
way we find numerically the distributions P(s, t) and Q(s, t) at any t for infinite N. Solving 10 6 
equations gives the evolution of these distributions for l<s<10 6 and S(t) = 1 — J2l=i P(s, t). 
The log-log plot (Fig. 2A) shows that the obtained S(t) dependence is well described by the 
power law: S oc 5 13 with S = t - t c , t c = 0.923207508(2), and small (5 = 0.0555(1) (which 
is very close if not equal to 1/18 = 0.05555...). Here we find t c as the point at which P(s) 
is power-law over the full range of s, see below. To check the correctness and precision of 
our calculations, we repeated them for ordinary percolation and obtained the classical results 
with the same presicion as for our model. Although the small exponent /3 makes it difficult to 
approach the narrow region of small S, fitting this data by the law + bd 13 , we find that is 
smaller than 0.005. This supports our hypothesis that the transition is continuous, but still does 
not prove it. Moreover, both our extensive simulations and numerical solution results clearly 
demonstrate that the analysis of the S(t) data cannot validate or rule out a discontinuity. 
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Figure 2B shows the evolution of the distribution P(s, t), which we compare with the evo- 
lution of this distribution for ordinary percolation, Fig. 2C. The difference is strong at t < t c , 
where the distribution for explosive percolation has a bump, but above t c the behaviors are sim- 
ilar. The distribution function Q(s, t) evolves similarly to P(s, t) in the full time range. At the 
critical point, we find power-law P(s) ~ s 1_r and Q(s) ~ s 3 ~ 2r in the full range of s (six orders 
of magnitude), where r = 2.04762(2), which is close to 2, as in (7,8,10), in contrast to r = 5/2 
for classical percolation. The first moments of these distributions, (s) P = J2 S sP(s) (the mean 
size of a finite cluster to which a randomly chosen node belong) and (s) Q = J2 S sQ(s), demon- 
strate power-law critical singularities (s) P ~ |5|~ 7p and (s) Q ~ |#|~ 7Q , where exponents 
7p = 1.0111(1) and 7 Q = 1.0556(5) both below and above the transition. Note that 7 P > 1 
in contrast to ordinary percolation, where the mean-field value of exponent 7 is 1. Figure 2D 
shows the set of time dependencies of P(s, t) for fixed cluster sizes (the time dependencies of 
Q(s,t) are similar). These dependencies strongly differ from those for ordinary percolation 
(Fig. 2E) in the following respect. The peaks in Fig. 2D for explosive percolation are below t c , 
while the peaks in Fig. 2E for ordinary percolation are symmetrical with respect to the critical 
point at large s. 

The inspection of these numerical results in the critical region (t < t c ) reveals a scal- 
ing behavior typical for continuous phase transitions, P(s,t) = s 1 ~ r /(s5 1 / cr ) and Q(s,t) = 
s 3 ~ 2t g(s5 1 / (7 ), respectively, where f(x) and g(x) are scaling functions, and a = (r — 2)//3, 
which is a standard scaling relation. One should stress that these functions are quite unusual. 
Figure 3 shows the resulting scaling functions and, for comparison, the scaling function for ordi- 
nary percolation. Remarkably, f(x) and, especially, g(x) are well fitted by Gaussian functions. 
These functions differ dramatically from the monotonously decaying exact scaling function 
e~ x I \/2x for ordinary percolation. Effective elimination of the smallest clusters in this merging 
process results in the minima of the scaling functions at x—0. On the other hand, the stunted 



6 



emergence of large clusters results in the particularly rapid decay of these functions at x 3> 1. 

The key point of our report is the following strict analytical derivation. We start from the 
fact that in this problem the cluster size distributions are power-law at the percolation threshold. 
We observed these power-laws over 6 orders of magnitude, and they were observed in works 
(7,8,10) for explosive percolation though in less wide range of s. Now we strictly show that 
if the cluster size distribution is power-law at the critical point, then this phase transition is 
continuous. We use the fact that in the critical region above the percolation threshold, the 
distributions Q(s) and P(s) at large s are proportional to each other, namely Q(s) = 2SP(s), 
see Eq. (1). This relation crucially simplifies our problem, since the resulting evolution equation 
for the asymptotics of the distribution in this region contains only P(s) and the relative size 
S(t) of the percolation cluster. As a result, Eq. (2) becomes very similar to that for ordinary 
percolation (the only difference is the presence of S(t) terms on the right-hand side), and so it 
can be easily analysed explicitly in the same way as for ordinary percolation (11). In this way, 
using a power-law critical distribution P(s) ~ s 1_r as an initial condition for this equation at 
t = t c , we find that the behavior of the distribution P(s) and S(t) in explosive percolation above 
the percolation threshold is qualitatively similar to that for ordinary percolation. Specifically, we 
show that the percolation cluster emerges continuously, and S oc 5 13 , where t = 2 + (5/(1 + 3/3) 
and so a = 1/(1 + 3/3). The obtained numerical values of exponents r, /3, and a agree with 
these two scaling relations and thus also confirm the continuous transition. So the results of this 
report are self-consistent. Our results are summarized in Table 1. Assuming a scaling form for 
the distributions gives 7 P = 1 + 2/3 and 7 Q = 1 + f3, which agree with our numerical solution of 
Eq. (2). Furthermore, applying standard scaling relations (1), we calculate the fractal dimension 
for this model, df = 2/a = 2(1 + 3/3), and the upper critical dimension, d u = df + 2/3 = 
2(1 + 4/3) (see Table 1). The latter determines the finite size effect: f c (oo) - t c (N) oc N~ 2/du , 
where 2/d u = 0.818(1). Interestingly, the obtained fractal and upper critical dimensions for 
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explosive percolation are less than 3. They are much smaller than those for ordinary percolation, 
which are 4 and 6, respectively. Our model is infinite-dimensional, that is above the upper 
critical dimension, where mean-field theories must work, which makes the observed smallness 
of exponent f3 particularly remarkable. We know no other model in statistical physics with a 
small /3 in such a situation. 

In the general case, each of the samplings in the process involves m > 1 nodes, and so 
the minimal cluster is selected from the m possibilities. In this case, the relation between 
between the distributions Q(s) and P(s) is substituted by a more general one, but the form of 
the evolution equation (2) does not change. We find that with increasing m, t c approaches 1 and 
j3 = t—2 ~ rrT l e~ m , that is, f3 rapidly decreases with m, but the transition remains continuous. 
Thus the critical exponents depend on m, and so the "explosive percolation" transitions have no 
universal critical behavior. 

We have shown that the "explosive percolation" transition is actually continuous. It is the 
smallness of the (5 exponent for the size of the percolation cluster that makes it virtually impos- 
sible to distinguish this phase transition from a discontinuous one even in very large systems. 
Indeed, suppose that N = 10 18 and (5 1/18. The addition of one link changes t by At = 1/N, 
which is the smallest time interval in the problem. Then a single step At = 1CT 18 from the per- 
colation threshold already gives S ~ (At)' 3 ~ 0.1. Other critical exponents and dimensions 
also differ radically from classical values. Furthermore, we have derived a complete set of 
scaling relations between the critical exponents for this problem, which were also supported 
by our numerical results. The real absence of explosion topples an already established view 
of explosive percolation. We believe, however, that, thanks to the observed unique properties 
of this phase transition, our findings make this new class of irreversible systems an even more 
appealing subject for further extensive exploration. 
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Table 1 : Threshold values, critical exponents, and fractal and upper critical dimensions for 
ordinary percolation and explosive (m = 2) one. Critical exponents for explosive percolation 
are expressed in terms of exponent /?: r = 1 + /3/(l + 3/3), a = 1/(1 + 3/3), 7 P = (3 — r)ja — 



1 + 2/3, 7q 


= (5 - 2t)/<t = 


1 + P,d f = 


2(1 + 3/3), d u = 2(1 + 


•4/3). 














r o" 


7p 


7q 


df 




Ordinary 


1/2 


1 


5/2 1/2 


1 




4 


6 


Explosive 


0.923207508(2) 


0.0555(1) 


2.04762(2) 0.857(3) 


1.111(1) 


1.0556(5) 


2.333(1) 


2.445(1) 
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Here we compare models of explosive percolation and present some details of our analytical 
calculations. 

1 Comparison between explosive percolation models 

Let us show that our model provides more efficient merging of small clusters than the model 
(3). So if the explosive percolation transition is continuous in our model, then the model (3) 
also has a continuous transition. 

Recall the selection rule in our model. At each step, choose two pairs of nodes, % and j, k 
and I, uniformly at random. Let they belong to clusters of sizes and Sj, s k , and s t , respectively. 
Connect the smallest cluster of the pair Sj and Sj, with the smallest cluster of the pair s k , and 
si. One can immediately see that this rule is equivalent to the following one. Let f(s, s') be 
an arbitrary monotonously growing function of its arguments, e.g., f(s, s') = ss'. Consider 
four possibilities to add a new link, ik, il, jk, and jl, and choose that one which provides 
the smallest value of f(s i: s k ), f(s u s{), f(sj, s k ), and f(sj, s t ), see Fig. 4A in the Supporting 
online material. 

In the model (3), the product selection rule is as follows. At each step choose nodes % j, k 
and / uniformly at random, and select from two possibilities, either connect % and k or connect 
j and /, choosing the pair with the smallest of the products SiS k and SjSi, see Fig. IB. 

Then the only difference between our model and the model (5) is that we select from four 
possibilities (comparison of SjS fc , SiS h SjS k and SjSi) while in the model (3) the selection is 
from only two possibilities (comparison of SiS k and SjSi). Therefore our choice guarantees 
merging of clusters with the product of sizes ss' which is equal or smaller than in the model 
(3). So our process should generate discontinuity even more efficiently than in the model (5). 
Consequently, if our model has a continuous phase transition, that the model (3) also must have 
a continuous transition. 
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2 Analysis of the evolution equation above the explosive per- 
colation transition 



Let us consider the master equation for our model 

^l^- = s J2 Q(u,t)Q(v,t)-2sQ(s,t), (3) 

°^ u+v=s 

where 

Q(s) = [Pcum(s) + PcuUs + 1) + 2S]P(s) = [2- 2P(1) - ... - 2P{s - 1) - P{s)]P{s). (4) 

(Here we only consider the case of m = 2.) We introduce generating functions for the distribu- 
tions P(s) and Q(s), 

oo 

p(z)=J2P(s)z s (5) 

and 

oo 

a(z) = J2Q(s)z s . (6) 

Note that p(z = 1) = 1 — Sand 17(2: = 1) = 1— S 2 , where Sis the relative size of the percolation 
cluster. Using these generating functions for the distributions we represent the master equation 
© in the following form: 

In ordinary percolation, cr(z, t) in this equation is substituted by p(z, t), namely 

In the critical region above the percolation threshold, Eq. © gives 

Q{s) 2SP(s) (9) 



14 



asymptotically at large s, which leads to a simple relation between the generating functions 
a(z) and p(z) in the range z close to 1. Indeed, 

1 _ s 2 - a(z) = £Q(s)[l ~ A =T, 2 SP(s)[l - z s \ = 2S[1 -S- p(z)}, (10) 

s s 

so 

1 - a{z) = 2S[1 - p(z) - S/2] (11) 
if z is close to 1 in the critical region at t > t c . One can check this relation at z = 1, namely, 

1 - = S 2 = 2S[1 - p(l) - S/2] = 2S(S - S/2). (12) 

Therefore, in the critical region above the percolation threshold, at z close to 1, the master 
equation takes a convenient form, 

3eM = M -( t)W ,, t )-i + sw/ 2 i^M. (B) 

Note that this equation essentially differs from Eq. ([8]) for ordinary percolation because of the 
terms S(t) on the right-hand side. Nonetheless Eq. (TT3T ) can be analysed in the same way as for 
ordinary percolation. 

Our numerical solution of Eq. © showed convincingly that at the critical point, the distri- 
bution P(s, t c ) is power-law, namely, at large s, P(s, t c ) = f(0)s 1 ^ T . Here /(0) is the critical 
amplitude for this distribution. This is also the value of the scaling function for this distribution, 
f(x = 0), see below. Let us show that, if at the critical point P(s, t c ) oc s 1_T , then Eq. (fT3l) has 
a solution with 1 — p(z = 1, t) = S(t) oc (t — t c )P in the critical region, which just means that 
the transition is continuous. The existence of this solution can be demonstrated in the following 
way. Let us assume that P(s, t c ) = /(0)s 1_r at large s, and S = B(t — t c Y at small t — t c and 
check whether this assumption is correct or not. Here we assume that /(0) and the exponent r 
are known (numerical solution gave /(0) = 0.04618(2) and r = 2.04762(2)), while B and the 
exponent f3 are to be found. 
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We use the power-law asymptotics of the distribution P(s,t c ) = f(0)s 1 ~ T as the initial 
condition for Eq. (fT3l) . This corresponds to the following singularity of the generating function 
atz = 1: 



1 — p(z, t c ) = analytic terms — /(0)T(2 — r)(l — z 



,T-2 



(14) 



Introducing e = (t — t c ) 2l3+l and x = Inz, we transform Eq. (fl~3l) to the following form: 

de 1 + 2/5 V 2 jcfe ^ J 

To solve this equation, we use the hodograph transformation approach. We pass from p = 
p(x, e) to x = x(p, e), which leads to a simple linear partial differential equation for x(p, e) and 
enables us to find the general solution 



In z 



8B 2 



1 + 2/3 



B 



c/3/(l+2/3) 



2 1 + /?/(! + 2/3) 



(16) 



where the function F(p) is obtained from the initial condition (fl4l) . which gives the solution 



In z 



8B 2 
1+2/3 



1-p- 



B (t- t c 



2 l+/3/(l+2/3) 



( H )W_ [/(0) ]-l/(r-2)| r(2 _ T) |-l/M) [l _ p ]l/(T-2) i 

(17) 

We set ,2=1, and taking into account 1 — p(l) = S = B(t — t^ 13 and comparing resulting 
powers and coefficients in Eq. (fT7l) . we obtain a relation between the critical exponents, 



t = 2 + 



(18) 



1 + 3/3' 

and express the critical amplitude B for the relative size of the percolation cluster in terms of 
the critical amplitude /(0) for the distribution P(s, t c ), 

B= A {r Lj[j 6r) [/(0)] 1 /( 7 - 3 ^|r(2-r)| 1 /( 7 - 3T ). (19) 

3 — T 

The relation (fT8l precisely agrees with our numerical results, r = 2.04762(2) and (3 = 0.0555(1). 
Substituting /(0) = 0.04618(2) (our numerical result) into Eq. (0 gives B = 1.075, which 
agrees with the corresponding value 1.080 obtained by solving the muster equation numerically. 
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Furthermore, the power-law distribution at the critical point can be justified strictly by using 
an equation for scaling functions. In the normal phase (t < t c ), we derived an equation for 
the scaling functions in this problem. This is a nonlinear integral-differential eigenfunction 
equation, where eigenfunctions are the scaling functions for P(s, t) and Q(s, t), see Eq. (|2Q|) . 
and a critical exponent, e.g., r, plays the role of an eigenvalue. This equation can be solved 
numerically, which is, however, a difficult task. We verified that the scaling functions and r, 
which we found numerically by solving Eq. ©, satisfy this equation. This shows that the 
distributions are power-law at the critical point. 

3 Relations between critical exponents 

Here we present a list of relations for critical exponents for explosive percolation (m = 2). 
In the critical region, the distributions P(s, t) and Q(s, t) have a scaling form: 

P(s,t) = s 1 - T f(s5 1 /°) 

Q(s,t) = s 3 ~ 2T g(s5 1 /°), (20) 

where 5 = \t — t c \. Note that the critical exponents below and above the transition are equal, 
while the scaling functions below t c differ dramatically from those above the transition. The 
exponent j3 of the size of the percolation cluster is expressed in terms of r and a as follows: 

0=— ■ (21) 
o 

For the first moments of the distributions P(s) and Q(s), critical exponents are 7 P and j Q , 
respectively. Namely, (s) P oc |<5|~ 7P and (s) Q oc |<5|~ 7P . These exponents are expressed in 



therms of r and a as follows: 



7p = (22) 

7. = (23) 
a 
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Note the relation between 7 P and 7 Q : 

lp + 1 = 2 7q . (24) 

Finally, we give the full set of critical exponents, the fractal dimension df = 2/ a and the 
upper critical dimension d u = df + 2(3 in terms of the exponent (3 in the case of m = 2: 

^ = 2 + IW (25) 
" = TW (26) 

7p = 1 + 2/3, (27) 

7q = 1 + /3, (28) 

d f = 2(1 + 3/3), (29) 

4 = 2(1 + 4/3). (30) 

4 Relation between the critical time t c and exponent r 

We can obtain approximate relations between t c and r or between the critical amplitude /(0) 
and r by applying the sum rule J2 S P( s ) = 1 at the critical point. Two estimates are possible. 
One can estimate P(s, t c ) by its asymptotics /(0)s 1_r , which gives 

/(0)C(r-l)«l, (31) 

where ((x) = Y%Li s ~ x ls m e Riemann zeta function. If r is close to 2, then ((r — 1) = 
l/(r — 2), so we have r — 2 w /(0). This estimate shows that the small values of r — 2 and /(0) 
are interrelated. Recall that we obtained numerically r = 2.04762(2) and /(0) = 0.04618(2). 

In the second estimate we use the following approximation: P(s,t c ) w P(l, t c )s 1-r . We 
find P(l, t) explicitly by solving the master equation ©, which gives 

p (M) = Tlb*. (32) 
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so we have 

n^C(r-i)«l. (33) 

Using t c = 0.923207508(2) obtained numerically, we find approximately r — 2 ps 0.05, that is, 
the exponent r is close to 2 when i c is close to I. 
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Figure 1 : Percolation processes. (A) In classical percolation, at each step, two randomly chosen 
nodes are connected by a new link. If these nodes belong to different clusters, these clusters 
merge. (B) In the "explosive percolation" model, at each step, two pairs of nodes are chosen at 
random, and for each of the pairs, the node belonging to the minimal cluster is chosen. These 
two nodes (and so their clusters) are connected by a new link. (C) The relative size of the 
percolation cluster S versus t (the ration of the number of links L and nodes N) obtained from 
the simulation of our model with N = 2 x 10 9 nodes (1000 runs). The staring configuration is 
N bare nodes. For comparison, the corresponding result for classical percolation is shown. (D) 
Log-log plot S versus t — t c . The data S of our simulation (black curve) is better fitted by the 
a(t — tc) 13 law (dashed blue curve) than by the So + b(t — law with S = 0.31. In the latter 
case, In (S — S ) (red curve) does not follow a linear law (dash-dotted line). 
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Figure 2: Numerical solutions of the evolution equations (2) for the infinite system. (A) Log- 
log plot S vs. t — t c . The slope of the dashed line is 0.0555. (B) The evolution of the distribution 
P(s) below (black lines) and above (red dashed lines) the percolation threshold for explosive 
percolation. The distribution at the critical point is show by the blue line. (C) The evolution 
of P(s) for ordinary (classical) percolation. (D) Dependence of P(s, t) on t for a set of cluster 
sizes s for explosive percolation. Numbers on curves indicate s. (E) P(s,t) versus for normal 
percolation. The insets show the P(s, t) curves for large values of s. 
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Figure 3: Scaling functions f(x) and g(x) for explosive percolation (t < t c ) and, for compar- 
ison, the exact scaling function fcp(x) = e~ 2x /y/2n for ordinary (classical) percolation. The 
blue dashed lines show the Gaussian fittings of the explosive percolation scaling functions. 




Figure 4: Comparison of the linking rules in our model (A) and in the model (3) (B). 
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